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Abstract. Using Regge calculus, we construct a Regge differential equation for 
the time evolution of the scale factor a(t) in the Einstein-de Sitter cosmology 
model (EdS). We propose two modifications to the Regge calculus approach: 1) we 
allow the graphical links on spatial hypersurfaces to be large, as in direct particle 
interaction when the interacting particles reside in different galaxies, and 2) we assume 
luminosity distance Dl is related to graphical proper distance D p by the equation 

Dl = (1 + z)\Jd^, ■ T)^, where the inner product can differ from its usual trivial form. 
The modified Regge calculus model (MORC), EdS and ACDM are compared using 
the data from the Union2 Compilation, i.e., distance moduli and redshifts for type 

la supernovae. We find that a best fit line through log [ — —J versus \ogz gives a 

V G P C / 

correlation of 0.9955 and a sum of squares error (SSE) of 1.95. By comparison, the 
best fit ACDM gives SSE = 1.79 using H a = 69.2 km/s/Mpc, Sl M = 0.29 and 0\ = 
0.71. The best fit EdS gives SSE = 2.68 using H a = 60.9 km/s/Mpc. The best fit 
MORC gives SSE = 1.77 and H a = 73.9 km/s/Mpc using R = A' 1 = 8.38 Gey and 
rn = 1.71 x 10 52 kg, where R is the current graphical proper distance between nodes, 
A^ 1 is the scaling factor from our non-trival inner product, and m is the nodal mass. 
Thus, MORC improves EdS as well as ACDM in accounting for distance moduli and 
redshifts for type la supernovae without having to invoke accelerated expansion, i.e., 
there is no dark energy and the universe is always decelerating. 
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1. Introduction 



The problem of cosmological "dark energy" is by now well known[lJ|2j[3j[4j[5j[6j. 
Essentially, redshifts and distance moduli for type la supernovae indicate the universe 
is in a state of accelerated expansion when analyzed using general relativistic 
cosmology[7][8][9]. Specifically, the distance moduli increase with increasing redshift 
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faster than predicted by general relativistic cosmology using matter alone. Until this 
discovery in 1998, the so-called "standard model of cosmology" was general relativistic 
cosmology with a perfect fluid stress-energy tensor and an early period of inflation. 
Since this leads to a decelerating expansion (except during the short, early inflationary 
period), something 'exotic' seemed to be required to account for the unexpectedly large 
distance moduli at larger redshifts, viz., dark energy that causes the universe to change 
from deceleration to acceleration at about z = 0.752 [9J. The new "standard model of 
cosmology," i.e., that with the most robust fit to all observational data (ACDM), simply 
adds a cosmological constant A to the Einstein-de Sitter cosmology model (Q M +Q^ = 1) 
and A then provides the mechanism for accelerated expansion, i.e., it provides the 
dark energy. The "problem" is that our best theories of quantum physics tell us the 
cosmological constant should be exactly zero[10j or something hideously large and 
neither of these two cases holds in ACDM. Thus, one of the most pressing problems 
in cosmology today is to account for the unexpectedly large distance moduli at larger 
redshifts observed for type la super novae [6]. 

The most popular attempts to explain the apparent accelerating expansion of 
the universe include quintessence [12] [T3] and inhomogeneous spacetime[lj[2j[3j[4j[14j 
(there are even combinations of the two[15j [16J). Although these solutions have 
their critics [TT], they are certainly promising approaches. Another popular attempt 
is the modification of general relativity (GR). These approaches, such as f(R) 
gravitvp] [15] [20] [2T] [22] [23] . have stimulated much debate [21] [25] [26], which is a 
healthy situation in science. Herein, we propose a new approach to the modification of 
GR via its graphical counterpart, Regge calculus. 

Specifically, we construct a Regge differential equation for the time evolution of the 
scale factor a(t) in the Einstein-de Sitter cosmology model (EdS), then we propose 
two modifications, both motivated by our work on foundational issues[27j[28j[29j. 
First, we allow spatial links of the Regge graph to be large (as defined below) in 
accord with 1) our form of direct particle interaction between sources in different 
galaxies and 2) the assumption that Regge calculus is fundamental while GR is 
the continuous approximation thereto. Of course, direct particle interaction in its 
original form would require a modification to general relativistic cosmology in and of 
itself |30j [3 lj[32j[33j[34j[35j. We are not concerned with saving direct particle interaction 
in its original form and, indeed, one needn't accept our version thereof to consider the 
modifications of GR proposed herein, i.e., empirical motivations suffice. Second, we do 
not assume that luminosity distance Dl is trivially related to graphical proper distance 
D p between photon receiver and emitter as it is in EdS, i.e., in EdS Dl = (1 + z)d p 
where d p is proper distance between photon receiver and emitter. There are two reasons 
we do not make this assumption. First, in our view, space, time and sources are co- 
constructed, yet Dp is found without taking into account EM sources responsible for Dl- 
That is to say, in Regge EdS (as in EdS) we assume that pressureless dust dominates 
the stress-energy tensor and is exclusively responsible for the graphical notion of spatial 
distance D p . However, even though the EM contribution to the stress-energy tensor is 
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negligible, EM sources are being used to measure the spatial distance D L . Second, in the 
continuous, GR view of photon exchange, one considers light rays (or wave fronts) in an 
expanding space to find Dl = (1 + z)d p . In our view, there are no "photon paths being 
stretched by expanding space," so we cannot simply assume D L = (1 + z)D p as in EdS. 
Indeed, we find the trivial EdS relationship between luminosity distance and proper 
distance holds only when D p is small on cosmological scales. In order to generate a 
relationship between D L and D p , we turned to the self-consistency equation KQ = J in 
our foundational approach to physics[28], where K is the differential operator, Q is the 
'fieldfj] and J is the source. Since we want a relationship between Dl and D pi the 'field' of 
interest is a metric h a p relating the graphical proper distance D p , obtained theoretically 
using no EM sources, to the luminosity distance Dl, obtained observationally via EM 
sources. The region in question (inter-nodal region between emitter and receiver) has 
metric r] a p given by ds 2 = —c 2 dt 2 + dD 2 , so the inner product of interest can be written 
Va/3 + ha/3 where the spatial coordinate is D p and h a p is diagonal. Since each EM source 
proper is not "stretched out" by the expansion of space, the spatiotemporal relationship 
between emitter and receiver is modeled per this inter-nodal region alone. Thus, unlike 
EdS, we have no a priori basis in our form of direct particle interaction to relate Dl 

to D pi so we begin with the assumption Dl = (1 + z)\Jl?' p ■ ~T? p = (1 + z)D p -J\ + h n , 

where D p = (0, D p ). 

The specific form of KQ = J that we used was borrowed from linearized gravity in 
the harmonic gauge, i.e., d 2 h a/ 3 = —16irG(T al3 — |^ Qj gT). We emphasize that h a p here 
corrects the graphical inner product 7] a p in the inter-nodal region between the worldlines 
of photon emitter and receiver, where 7] a p is obtained via a matter-only stress-energy 
tensor. Since the EM sources are negligible in the matter-dominated solution, we have 
d 2 h al3 = to be solved for h\\. Obviously, h u = is the solution that gives the trivial 
relationship, but allowing hu to be a function of D p allows for the possibility that Dl 
and D p are not trivially related. We have hu = AD P + B where A and B are constants 
and, if the inner product is to reduce to r\ a p for small D p , we have B = 0. Presumably, 
A should follow from the corresponding theory of quantum gravity, so an experimental 
determination of its value provides a guide to quantum gravity per our view of classical 
gravity. As we will show, our best fit to the Union2 Compilation data gives A^ 1 = 
8.38 Gey, so the correction to r/n is negligible except at cosmological distances, as 
expected. Essentially, we're saying the dark energy phenomenon is an indication that 
A 7^ so that one cannot simply assume the distance Dl measured using EM sources 
corresponds trivially to the graphical proper distance D p even though the EM sources 
contribute negligibly to the stress-energy tensor. 

One might also ask about distance corrections per h 0Q , i.e., as regards redshift, 
but since redshift distances are fractions of a meter one wouldn't expect h 0Q to be of 
consequence here. Of course, there is the issue of origin of redshift in our approach, since 

| The interested reader is referred to section 3 of reference |28j for an explanation of how our notion 
of a "held" is consistent with our notion of direct particle interaction. 
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typically cosmological redshift is understood to occur between emission and reception [ 
while clearly it must occur during emission and reception in our view. While we don't 
have photons propagating through otherwise empty space between emitter and receiver, 
we do relate the reception and emission events in null fashion through the simplices 
spanning the inter-nodal region between emitter and receiver. Using the metric in each 
simplex ds 2 = —c 2 dt 2 + dD 2 , as above, we have dD p = adx, just as in EdS, although t is 
not proper time for the nodal observers as it is in EdS. This difference in t is accounted 
for in the computation of D p where it has a small effect for the range of data in the 
Union2 Compilation^. Likewise, we do not find that it leads to a significant difference 
in scale factor at time of emission a e as a function of z for the data range in question. 
Not surprisingly, when we compute the redshift graphically we find it is equivalent to 

the special relativity (SR) result, i.e., z + 1 = \ -h— - Y e Jf^})~ Y^fl where V e is the 

V ^ V (1-V e /C)(l-V r / C ) 

velocity of the emitter at time of emission in the (l+l)-dimensional inter-nodal frame 
and V r is the velocity of the receiver at time of reception. Using this form of redshift 
in the EdS model and comparing the result to the use of cosmological redshift in EdS, 
we find there is no significant difference between the two results for distance modulus 
fi versus redshift z well beyond the range of the Union2 Compilation (z < 2, see Figure 

|2l). Therefore, we use cosmological redshift a e = for the computation of D„, since 

1 + z 

cosmological redshift is far simpler than the graphical alternative. 

While these modifications are motivated by our work on foundational issues, their 
specific mathematical instantiations are herein aimed at explaining dark energy. Since 
this is our first foray into modified Regge calculus (MORC), the specific approaches 
required for explaining other GR phenomena, e.g., the perihelion shift of Mercury, 
remain to be seen. A defense of MORC will not be undertaken here, interested readers 
are referred to our earlier work cited above, but a couple comments are perhaps in 
order. First, the graphical lattice used herein to obtain a(t) clearly violates isotropy 
and is not to be understood as a literal picture of the distribution of matter in the 
universe, e.g., galactic clusters, voids, etc. In a sense, the graphical lattice we use is 
no coarser an approximation than the continuum counterpart it is designed to replace, 
i.e., the featureless perfect fluid model of EdS where there is absolutely no structure. 
Rather, the graphical lattice simply provides a 'mean' evolution for the scale factor 
a(t) in the equation for D p . Second, the goal of such idealized models is to attempt to 
isolate 'average' geometric and/or material features of cosmology which broadly capture 
kinematic properties of the universe as a whole. Only when such models show some 
initial success are explorations into departures from their simplistic structure motivated, 
e.g., the inhomogeneous spacetime models cited above. Thus, the model we present 
herein was designed merely to test the possibility of replacing the continuous EdS 
cosmology with a discrete, graphical counterpart based on our form of direct particle 

cdt 

§ There is another difference between d p and D p as computed using d\ = — that must be considered. 
This wiil be explained in section [2] 
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interaction (again, for reasons unrelated to dark energy). Only upon some success of 
this initial test, i.e., improving the EdS fit to the type la supernova data, should we 
proceed to address the commensurate questions and implications of this approach (as 
outlined briefly in section H] of this paper). We believe the results presented herein 
establish precisely "some initial success" and therefore justify further exploration into 
this idea. 

We begin in section [2] with an overview of Regge calculus and present our temporally 
continuous, spatially discrete Regge EdS equation for the time evolution of the scale 
factor a(t) and the commensurate equation for proper distance between photon emitter 
and receiver D v in a direct inter-nodal exchange. As we will see, the spatially discrete 
Regge EdS equation for the time evolution of the scale factor a(t) reproduces that of EdS 
when spatial links are small. Spatial links are "small" when the 'Newtonian' graphical 
velocity v between spatially adjacent nodes on the Regge graph is small compared to 

c, i.e., f— j < 1. In that case the dynamics between adjacent spatial nodes is just 
Newtonian and the evolution of a(t) in Regge EdS is equal to that in EdS. Deviations 
in the evolution of a(t) between Regge EdS and MORC turn out to be small (see 
Figure ED- Thus, the modification of Regge evolution plays a relatively minor role 
in the MORC fits. Rather, as we will show, the major factor in improving EdS is 



Dl = (1 + z)D p — > Dl = (1 + z) A/ D p ■ D p . Since Regge EdS should give EdS when 
used as originally intended [37], the proposed mechanism for EM coupling in MORC 
differs from that in Regge calculus. When v ~ 2c Regge EdS encounters the "stop 
point" problem[38j[39j[40j, i.e., the backward time evolution of a(t) halts, so a(t) has a 
minimum and there is a maximum value of z for which one can find D p . Of course, this 
is not a real problem for Regge EdS if one is simply using it to model EdS, since one 
can regularly check v in the computational algorithm and refine the size of the lattice 
to ensure v remains small. However, in our case the graphical approach is fundamental, 
so lattice refinements are not mere mathematical adjustments, but would constitute 
new 'mean' configurations of matter. Of course, such refinements are certainly required 
in earlier cosmological eras, but one would expect there exists a smallest spatial scale 
(associated with a smallest nodal mass) so that eventually (evolving backwards in time) 
v ~ 2c could not be avoided and the minimum a(t) would be reached. Thus, there are 
significant deviations from our use of Regge calculus and its (originally intended) use as 
a graphical approximation to GR. 

In section |3] we present the fits for EdS, MORC, and ACDM to the Union2 
Compilation data, i.e., distance moduli and redshifts for type la supernovaejH] (see 
Figures H] and [5]) . We find that MORC improves EdS as much as ACDM in accounting 
for distance moduli and redshifts for type la supernovae even though the MORC universe 
contains no dark energy is therefore always decelerating. While we do not need to invoke 
dark energy, we do propose modifications to classical gravity. Thus, it is a matter of 
debate as to which approach (ACDM or MORC) is better. 

Of course, the success of MORC in this context does not commit one to our 
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Figure 1. (a) Tessellated sphere and (b) two "flattened" trapezoids (green) from the 
sphere. 



foundational motives. In fact, one may certainly dismiss our form of direct particle 
interaction and simply suppose that the metric established by EM sources deviates 
from that of pressureless dust at cosmological distances in a graphical approach to 
gravity. Since motives are not germane to physics, we will not present arguments for our 
foundational motives here. Abandoning our motives but keeping the MORC formalism 
would simply result in a situation similar to that in ACDM where a cosmological constant 
is added to EdS for empirical reasons. That is, one could simply view MORC as 
a modification of Regge calculus for empirical reasons without buying into our story 
about direct particle interaction and co-constructed space, time and sources. Motives 
notwithstanding, we believe our MORC formalism may provide creative new approaches 
to other long-standing problems, e.g., quantum gravity, unification, and dark matter. 
We conclude in section H] by briefly outlining future directions and challenges for this 
research program. 



2. Overview of Regge Calculus 

Regge calculus is typically viewed as a discrete approximation to GR where the discrete 
counterpart to Einstein's equations is obtained from the least action principal on a 4D 
graph[37J[42j|43j[44j. This generates a rule for constructing a discrete approximation 
to the spacetime manifold of GR using small, contiguous 4D Minkowskian graphical 
'tetrahedra' called "simplices." The smaller the legs of the simplices, the better 
one may approximate a different iable manifold via a lattice spacetime (Figure [1]). 
Although the lattice geometry is typically viewed as an approximation to the continuous 
spacetime manifold, it could be that discrete spacetime is fundamental while "the usual 
continuum theory is very likely only an approximation [45] " and that is what we assume. 
Curvature in Regge calculus is represented by "deficit angles" (Figure [I]) about any 
plane orthogonal to a "hinge" (triangular side of a tetrahedron, which is a side of a 
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Figure 2. Comparison of cosmological redshift (blue) and graphical special relativistic 
redshift (red) using EdS. The two curves begin to be resolved at z = 6. 



simplest]]]) , so curvature is said to reside on the hinges. A hinge is two dimensions less 
than the lattice dimension, so in 2D a hinge is a zero-dimensional point (Figure [Q. 

The Hilbert action for a vacuum lattice is Ir = — } EiAi where <7j is a triangular 



hinge in the lattice L, is the area of ex, and £j is the deficit angle associated with e^. 

SIr 



The counterpart to Einstein's equations is then obtained by demanding 



where 



£j is the squared length of the j lattice edge, i.e., the metric. To obtain equations 
in the presence of matter-energy, one simply adds the matter-energy action Im to Ir 

and carries out the variation as before to obtain — ^ = — — ^ ESI - The LHS becomes 



51 



R 



51) 



1 

167T 



Si cot Qij where 0^ is the angle opposite edge £j in hinge One 



finds the stress-energy tensor is associated with lattice edges, just as the metric, and 
Regge's equations are to be satisfied for any particular choice of the two tensors on the 
lattice. The extent to which Regge calculus reproduces GR has been studied[47j[48j[49j 
and general methods for obtaining Regge equations have been produced[50], but these 
results are of no immediate concern to us because we simply seek the Regge counterpart 
to a specific GR equation, i.e., a Regge differential equation for the time evolution of 
the scale factor a(t) in EdS. Whether or not we obtain said equation will be clear by 
virtue of its ability to track the analytic EdS solution in the proper regime, so we will 
not have to delve into issues associated with the 'accuracy' of Regge calculus in general. 



Our hinges are triangles, but one may use other 2D polyhedra. 



Modified Regge calculus as an explanation of dark energy 
2.1. Regge EdS Equation and MORC 



Following Brewin[5T] and Gentle[52], we take the stress energy associated with the 
worldlines of our particles to be of the form 



YIGm 
c 2 (icAt) 



so our Regge equation is 



\2iR(a n + a n+ i) 




2 cos 



-i 



V3(f) 2 ( 


'a n+1 -a n - 


)* + 4 






' a n+l~ a n * 




. At , 



\2iGm 



c 3 At 

Multiplying both sides of (jTJ) by —icAt/12 and letting v = R(a n+ i — a n )/At gives 



(i) 



i2(a„ + a n+1 ) = -3- (2) 



c 2 



If At — > 0, then v can be regarded as a 'Newtonian' velocity and i?(a n + a n+ i) can be 
replaced by 2r, where r is the graphical proper distance between two adjacent vertices 
on the lattice. Equation (J5J) then becomes 



n _ cos -i f v 2 /* 2 \ _ 2 cos- 1 I vg55 



(3) 



V^V^T4 2rc 2 

which we emphasize is unmodified Regge calculus. If v 2 /c 2 1, then a power series 
expansion of the LHS of Equation (jHJ) gives 

^ + P) 4 > (4) 
4c 2 \c) 2rc 2 1 J 

i) 2 Gm 

Thus, to leading order, our Regge EdS is EdS, i.e., — = , which is just a Newtonian 

2 r 

conservation of energy expression for a unit mass moving at escape velocity v at distance 
r from mass m. To better understand the relationship between Regge EdS and EdS, 
we note that in EdS any comoving observer A can ask, "What is the proper time rate 
of change of proper distance for comoving observer B at a proper distance r away from 

i) 2 Gm 

me today?" The answer is precisely v given by the EdS equation — = , where 

2 r 

m is the mass contained inside the sphere of radius r centered on observer A. In EdS 
the matter is distributed uniformly throughout space so the mass m inside sphere of 
radius r goes as r 3 , thus v oc r on spatial hypersurfaces in the EdS equation, so there 
is no limit to how large v is in this expression, it's Newtonian. In Regge EdS, v is the 
relative 'Newtonian' velocity of spatially adjacent nodes of mass m. In our view, photon 
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exchanges occur in direct node-to-node fashion, but solving for a Regge graph between 
all galaxies in the universe is of course unreasonable. Instead, we use Equation (J3J) to 
provide a 'mean' a(t) for the computation of graphical proper distance D p between any 
two photon exchangers, as in EdS, i.e., 

f to dt f 1 da 

proper distance = \ e = c / — = c / — (5) 

he a Ja e a a 

We then compute D p as a function of z by using Equation (j3J) obtained from the 

'mean' graph. However, before we continue there are two issues that we need to address 

regarding Equation flSJ). 

First, while it is true that cdt = adx for a null path in a simplex and the null path 

cdt 

will cross all values of x between emitter and receiver, the sum of d\ = — will not 

a 

equal Xe, i-e., the radial coordinate of the emitter. That's because the lines of constant 

X are tilted in the simplices (Figure EJ), so there is a fraction of d\ (given by A in 

cdt 

Figure EJ) that is not accounted for by — . This A is positive on the emitter's side of 

a 

the simplex and negative on the receiver's side, but the A sum on the two sides won't 
cancel out exactly, since the extent of constant-x tilt is reduced during the expansion. 
The correct equation for the graph is 

— r( i+ T(f-8)i 

where V is the SR velocity of the emitter or receiver as a function of time and relates 
to our 'Newtonian' v per 

c y/1 + v 2 /4c 2 

To simplify the analysis and obtain an estimate of how much A contributes, we use 
EdS with z = 2 and H a = 70 km/s/Mpc. From EdS we have a(t) = 



fX 2/3 



— = 1 - (t 1/3 - tl /3 ), and — = . For z = 2 and H = 70 km/s/Mpc 

Xe Xe V ' 2c 3ct 2 J V/3 

we have t Q = 9.31 Gy, t e = 1.79 Gy, and Xe = 11.81 Gey Using these values in Equation 
(jUJ) we find (iteratively) Xe = 12.189 Gey This increases /i (Equation f|T3|) below) by 
0.069 at z = 2 where /i is slightly greater than 44 (Figure [5]). This increase adds 0.0137 

to log f — —J in our curve fitting, which amounts to a 1.3 percent increase at z = 2. 
VGpcy 

This change is only 0.75 percent at z — 0.5 and 0.004 percent at z = 0.1. Thus, given 
the scatter in the data, we will ignore this correction. 

Second, in EdS, the scaling factor at emission is related to the redshift by 

a P = . In EdS, this redshift is understood to occur while the radiation is in transit 

l + z 

between emitter and receiver [36]. This "cosmological" redshift can be understood in 
the graphical picture to result from the fact that dt in EdS runs along lines of constant 
X and these lines are tilted away from the center of the simplex towards its nodal 
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dt 




Figure 3. Lines of constant x are tilted away from midpoint of simplex towards 
emitter and receiver. 



f ° cdt 

worldlines as discussed above (Figure |3]). That is, A = in EdS so \ e = / — holds 

Jt e « 

exactly. Thus, two EdS null paths eminating from different points on a spatial link 
have their proper distance of separation increase from simplex to simplex. However, 
as explained above, our dt is perpendicular to the spatial links so the null paths of 
successive emissions do not increase proper distance separation when traced through 
the simplices, i.e., redshift occurs entirely at emission and reception. Thus, relating 
successive events along the emitter's worldline in null fashion to events on the receiver's 
worldline, it is not surprising that we find the time delay between successive reception 
events as related to the temporal spacing of the emission events is that given by SR, 
i.e., 



z + 1 



V e /c)(l + V r /c) 




V e /c)(l-V r /c) (8) 

where V e is the SR velocity of the emitter at time of emission in the (l+l)-dimensional 
inter-nodal frame and V r is the SR velocity of the receiver at time of reception. Again, 
these SR velocities relate to our graphical 'Newtonian' v per Equation (J7J). As above, we 

simplify the analysis using the EdS equation for a(t) and find v r = XeH Q and v e = e ° 



where, again, \ e is the comoving coordinate of the emitter with the receiver at the 
origin. We need to find ^J~a~ e as a function of z, then substitute into the equation for 
proper distance between photon exchangers in EdS 



dp = H: 



(9) 



Even with the simplifications, the process gets messy and ultimately was solved 

XeH Q 



numerically. Since a D 
and we find 



1, we have d p = Xe (as assumed in Equation [5]). Let x 



(A + l] 



{A -I) 



2c 



(10) 
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where 



A 



x/- 



X 



vT + x 2 + x 



Thus, Equation is x = 1 



land gives 



(A — l) 2 

A 2 - 2A + 1 - 2x A 2 + AAx -2x + A 2 x 2 + x 2 - 6Ar 2 







(12) 



We then solve Equation f lT2|) numerically for x as a function of 2 and compare with 

the EdS version, i.e., x = 1 , to obtain Figure [2] where we see that there is 

no significant difference between the two results well beyond the range of the Union2 
Compilation (z < 2). 

Since these two differences between MORC and EdS do not result in any significant 

difference in our fit to the data of interest, we simply use Equation fl5]) with a e = -^—^ — 
to compute D p . However, there is one additional difference between d p and D p when 
using Equation fl5]) that we will not ignore. We will address this additional (simple) 
correction in the following section where we fit EdS, MORC, and ACDM to the Union2 
Compilation. 



3. Data Analysis 

The Union2 Compilation provides distance modulus fi and redshift z for each supernova. 
In order to find /i versus z for each model, we first find proper distance as a function of 
z, then compute the luminosity distance D L , and finally 



\i = 5 log 



lOpc 



(13) 



For EdS we have Equation for d p , so the only parameter in fitting EdS is H Q . For 
ACDM we have a 



H, 



M 



Vt^a 2 where Vt M + = 1. Plugging this into Equation 



we obtain 




cos 



yrt^ — (a/3 — 1) s/nx 

3 M^+ (v^+i) a/HI 

+ i) v / ix:-(v / 3-i)v / n: 



2W3 N 
4 

+ v / 3' 



4 



(14) 



f <p -1/2 

where F(cj)\m) = (l — msin 2 ^) d6 is the elliptic integral of the first kind. Thus 
Jo 

there are two fitting parameters for ACDM, H and either VLm or Vt\. For MORC, 
Equation ([3]) gives us a(d) rather than d(a), so we modify Equation (|5]) to read 



D n = R 



bf(b) 



b 2 

l + —db 
4 



(15) 
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where b = Ra/c, 

m = ^E±i — „ (16) 

2["-cos-'(lfe)-2cos-'(^±i 
and 61 and b e respectively solve 

1 = ~TE^ hl ) and ° e = 

Crit (Tit 



6 2 

The factor W 1 + — is the correction needed to adjust the time dt in Equation <^ to 

proper time dr of the nodal worldlines. [This is the "one additional difference between 
dp and D p when using Equation ([5])" alluded to at the end of section [2j] Equation (jSJ) is 
then solved numerically for D p and Dl = (1 + z)D v sJ\ + v4.D p as explained in section 
[TJ There are three fitting parameters for MORC, the inter-nodal coordinate R on the 
'mean' graph, the nodal mass m on the 'mean' graph, and A~ l from h\\. Specifying m 

and R is equivalent to specifying H Q in EdS, i.e., H Q = \j - _ - in EdS with p given 

4 o 

by the graphical values of R and m per -irR p = m. Thus compared to EdS, MORC 

3 

(as with ACDM) has one additional fitting parameter A -1 , which presumably will be 
accounted for ultimately by the corresponding theory of quantum gravity. 

As mentioned above, we fit these three models to the Union2 Compilation data 
(see Figures H] and [5]) . In order to establish a statistical reference, we first found that 

D L 



a best fit line through log — — versus logz gives a correlation of 0.9955 and a sum 

VGpcy 

of squares error (SSE) of 1.95. EdS cannot produce a better fit than this best fit line. 
The best fit EdS gives SSE = 2.68 using H a = 60.9 km/s/Mpc. A current (2011) "best 
estimate" for the Hubble constant is H a = (73.8 ± 2.4) km/s/Mpc [53]. Both MORC 
and ACDM produce better fits than the best fit line with better values for the Hubble 
constant than EdS. The best fit ACDM gives SSE = 1.79 using H = 69.2 km/s/Mpc, 
tt M = 0.29 and fi A = 0.71. This best fit ACDM is consistent with its fit to the WMAP 
data using the latest distance measurements from BAO and a recent value of the Hubble 
constant[54j. The best fit MORC (case 1, Table 1) gives SSE = 1.77 and H a = 73.9 
km/s/Mpc using R = A^ 1 = 8.38 Gey and m = 1.71 x 10 52 kg. Given the scatter in the 
data, MORC and ACDM produce essentially equivalent fits, clearly superior to EdS. 

The "stop point" value of z in the MORC best fit is only 2.05, so we expect the 
Regge evolution deviates discernibly from the EdS evolution in this trial. To check 
this, we compared the Regge model using the best fit parameters and h\\ = with 
its EdS counterpart. As explained above, the EdS counterpart to a Regge graphical 

result is obtained by using H Q = \j - _ - in EdS with p given by the graphical values 
4 

of R and m per -nR 3 p = m. The top graph in Figure [6] shows there is in fact a 
3 

discernible difference between the Regge and EdS evolutions, and the EdS value of H a 
obtained per R and m in this trial is 68.5 km/s/Mpc, which is significantly lower than 
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log(D L /Gpc) 




log(z) 



Figure 4. Plot of transformed Union2 data along with the best fits for linear regression 
(gray), EdS (green), ACDM (blue), and MORC (red). 



H = 73.9 km/s/Mpc found in MORC. In fact, the twenty trials with the lowest SSE 
values (cases 1-20, Table 1) have "stop point" z less than 10, so Regge evolution, as 
distinct from EdS evolution, does come into play. However, Regge evolution tracks EdS 
evolution when "stop point" z is as small 9.98 (see bottom graph in Figure E]) as is 
true in case 21 of Table 1. And, SSE = 1.78 for case 21 is still comparable to SSE 
= 1.79 of the best fit ACDM. The only casuality in the higher "stop point" z trials 
is H Q , which is lowered when Regge evolution tracks EdS evolution. However, the H Q 
= 71.2 km/s/Mpc in case 21 is still comparable to H Q = 69.2 km/s/Mpc for the best 
fit ACDM. Thus, the Regge evolution plays a relatively minor role in the MORC fits. 



Since we used the cosmological redshift, Xe — — , and the Regge evolution played 



a minor role in the MORC fits, we conclude that the major factor in improving EdS 



is Dl = (1 + z)D p — > Dl = (1 + z) y D p ■ D p . Again, given the scatter of the Union2 
Compilation data, we consider any of the 35 MORC results in Table [T], where SSE < 
1.78 and H ranges (69.9 75.3) km/s/Mpc, equivalent to the best fit ACDM. 

4. Discussion 

We have explored a modified Regge calculus (MORC) approach to Einstein-de Sitter 
cosmology (EdS), comparing the result with ACDM using the Union2 Compilation 
of type la supernova data. Our motivation for MORC comes from our approach to 
foundational physics that involves a form of direct particle interaction whereby sources, 
space and time are co-constructed per a self-consistency equation. Accordingly, since 
EM sources are used to measure luminosity distance Dl but are not used to compute 
graphical proper distance D p , we did not expect D p to correspond trivially to the 
luminosity distance Dl, i.e., we did not assume Dl = (1 + z)D p . Rather, we assumed 



a more general relationship Dl = (1 + z)y D p ■ D p where the inner product employed a 
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Figure 5. Plot of Union2 data along with the best fits for EdS (green), ACDM (blue), 
and MORC (red). The MORC curve is terminated at z = 1.4 in this figure so that the 
ACDM curve is visible. 

correction to the inter-nodal graphical metric, T) a p — > r] a p + h a @ with spatial coordinate 
D p and h a p diagonal, so that Dl = (1 + z)D p y/l + h n . The method used to find hn 
was a form of our self-consistency equation KQ = J borrowed from the homogeneous 
linearized gravity equation in the harmonic gauge, i.e., d 2 h a p = 0. While hn = is the 
solution typically used, we allowed hn to be a function of D p which gave hn = AD P + B 
where A and B are constants. Since we wanted the inner product to reduce to rj a p for 
small D p , we set B = 0. Our best fit MORC (case 1, Table 1) gave A' 1 = 8.38 Gey, so 
the correction to rjn is negligible except at cosmological distances, as expected. 

We found that in the context of the Union2 Compilation data MORC improved 
EdS as well as ACDM without having to employ dark energy. That is, the MORC 
universe evolves per pressureless dust and is always decelerating yet it accounts for 
distance moduli versus redshifts for type la supernovae as well as ACDM. Of course, 
this does not commit one to our foundational motives. In fact, one may certainly dismiss 
our form of direct particle interaction and simply suppose that the metric established 
by EM sources deviates from that of pressureless dust at cosmological distances; we 
did not present arguments for our foundational motives here. Abandoning our motives 
while keeping the MORC formalism would simply result in a situation similar to that 
in ACDM where a cosmological constant is added to EdS for empirical reasons, i.e., 
Regge calculus was modified to account for distance moduli versus redshifts in type la 
supernovae. Motives notwithstanding, MORC's empirical success in dealing with dark 
energy gives us reason to believe this formal approach to classical gravity may provide 
creative new techniques for solving other long-standing problems, e.g., quantum gravity, 
unification, and dark matter. 

In order to explore this possibility, we need to check MORC against the 
Schwarzschild solution, where experimental data is well established and GR is well 
supported. While tests of the Schwarzschild solution have been conducted on spatial 
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Figure 6. Top graph shows Regge evolution (red) without h\x correction and EdS 
evolution (green) for case 1 Table 1 where the "stop point" z is 2.05. The bottom 
graph makes the same comparison for case 21 Table 1 where the "stop point" z is 9.98. 

scales much smaller than the cosmo logical scales where we found a correction to EdS, 
it has been shown that the simplices must be small in order to reproduce the GR 
redshift and the perihelion precession of Mercury in the Schwarzschild solution|55j [56J. 
Thus, we need to verify that MORC is consistent with the Schwarzschild solution per 
observational data. We might refine our study of MORC cosmology, but we feel the 
easiest way to test MORC is via the Schwarzschild solution where perhaps the issue of 
dark matter can be addressed in a fashion similar to dark energy in EdS. If by chance 
we are able to construct a MORC for the Schwarzschild solution that passes empirical 
muster, we would then consider the more general issue of an action for modified Regge 
calculus in order to consider new approaches to quantum gravity and unification. Given 
the level of uncertainty involved in the next step alone, we won't speculate further. 
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Table 1. Table of 35 trials that produced the best fits for MORC. Column R is 
X in R = 10 x m. Column p is X in p = X x 10~ 27 kg/m . Column A^ 1 is X in 
A^ 1 = 10 x m. The other columns are self explanatory. 
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